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We clarify different definitions of the density matrix by proposing the use of different names, 
the full density matrix for a single-closed quantum system, the compressed density matrix for the 
averaged single molecule state from an ensemble of molecules, and the reduced density matrix 
for a part of an entangled quantum system, respectively. We show that ensembles with the same 
compressed density matrix can be physically distinguished by observing fluctuations of various 
observables. This is in contrast to a general belief that ensembles with the same compressed density 
matrix are identical. Explicit expression for the fluctuation of an observable in a specifled ensemble 
is given. We have discussed the nature of nuclear magnetic resonance quantum computing. We 
show that the conclusion that there is no quantum entanglement in the current nuclear magnetic 
resonance quantum computing experiment is based on the unjustifled belief that ensembles 
having the same compressed density matrix are identical physically. Related issues in quantum 
communication are also discussed. 
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1. INTRODUCTION 

1 . In a quantum system, the wave function specifies all the physical properties of the system. The wave function 
provides the maximal information. A system with a wave function description is said to be in a pure state. In 
contrast sometimes there is incomplete information about a quantum system and this may be more common. 
Von Neumann and Landau proposed to use the density matrix to describe the state of a quantum system in 



this situation 



It has been widely used nowadays, and is presented in standard textbooks and monographs 



waoays, and is 



on quantum mechanics, for instance in Refs. |3|, U, |5|, |a, l3| ■ With the development of quantum information 
in recent years, there have been renewed interests in the density matrix formalism. It is the basic tool for the 
study of quantum entanglement, channel capacity in quantum communication and so on. 

2. Though the density matrix has been proposed for nearly 80 years, there are still confusions and sometime 
disputes on matters related to it. Some are caused by a confusion in the terminology, but others arc true 
disputes of the fundamental nature. A notable example is the question of distinctness of ensembles having 
the same density matrix. These matters are becoming more acute with the progress in quantum computation 
research and to cause true concerns. Whether current liquid nuclear magnetic resonance (NMR) quantum 
computation is genuine quantum mechanical or just a classical simulation is just one such problem. 

3. In this article, we shall carefully inspect the definition of the density matrix, and that of the mixed state. This 
will remove a large number of disagreements. We propose to use three different terminologies for the density 
matrix in different situations: the full density matrix to describe the state of a single, closed quantum system, 
it always describes a pure state and is equivalent to the wave function; a reduced density matrix for the state 
of part of a coupled quantum system, it usually describes a mixed state; a compressed density matrix for an 
ensemble of independent quantum systems. After the classification of the density matrix, we will examine the 
distinctness of ensembles having the same compressed density matrix. We will show that ensembles having the 
same compressed density matrix can be distinguished physically. We show that this does not contradict one of 
the fundamental postulates of quantum mechanics that the density matrix completely specifies all the properties 
of a quantum ensemble, as stated by Asher Peres^. We then apply this result to the question of the nature of 
current NMR quantum computation, and show that the conclusion that there is no entanglement in the current 
NMR quantum computing is based on the unjustified belief that ensembles having the same compressed density 
matrix are physically identical, and hence is invalid. The implication of this result in quantum key distribution 
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and other areas of quantum information are also discussed. 

2. DENSITY MATRIX AND MIXED STATE 

4. First we introduce some notations in order to avoid confusion because different people use different definitions. 

A molecule in this article refers to a single quantum system such as a qubit, a composite quantum system and 
so on. 

An ensemble in this article is defined as a collection of N non-interacting molecules of the same kind in the 
same or different quantum states. An ensemble contains N non-interacting molecules, and there are Ni number 



of molecules in state i — ^, ■ ■ ■ , k, X^iLi = 



This definition of an ensemble has been used by Shankar, Merzbacher, d'Espagnat, Von Neumann 
A molecule is a quantum system which may contain a complex internal structure, for instance an atom, a 
molecule, a Bose-Einstein condensate all may be viewed as a molecule in an ensemble. The molecules are of 
the same kind, but they are independent of each other and can be distinguished, for example by their position, 
hence they are not identical particles. Usually they are far apart, so that their de Broglie waves rarely overlap. 

We note that the definition of ensemble by Asher Peres is different from ours In his definition, an ensemble 
is defined by a probability distribution and there are infinite number of molecules in an ensemble. 



It is important to perceive the differences in the definition of ensemble in statistical physics and quantum 



mechanics. In statistical mechanics 



I, an ensemble is a collection of infinite number of the system under 
identical physical constraints. For instance, a micro-canonical ensemble is specified by the particle number 
N and the total energy E. Then a system with Ni number of molecules in state \ipi), and iV iVi + iV2 -|- 
■ ■ ■ + Ni + ■ ■ ■ satisfying the given constraints on the particle number and total energy will be a member in the 
statistical ensemble. Whereas in quantum mechanics, the system of N molecules itself is now an ensemble, each 
molecule is treated as a system and the density matrix p describes the state of averaged molecule in this system. 
This different comprehension of ensembles in statistical mechanics and quantum mechanics causes many of the 
misunderstanding. 

A coupled quantum system is a quantum system with several constituents. For instance, a molecule is 
a coupled quantum system, and each atom in a multi-atom molecule is a constituent. The constituents in a 



coupled quantum system may be confined, for instance the atoms in a molecule, or far away apart, for instance 
two entangled photon pairs forming a coupled quantum system. 

By a sampling measurement of observable fl, we mean the following. We select randomly one molecule 
from an ensemble and makes an f2 measurement on this molecule. Then the ensemble is restored to its original 
state, namely Ni molecules in state |V-'i), and we randomly select a molecule again and makes an Q measurement 
on the molecule. We repeat this process sufficient number of times. Then the result of the sampling measurement 
of 17 is the averaged value of these individual measurements, and it is denoted by (il). In the density matrix 
formalism, the expectation value of an observable Q in an ensemble, (f2) = Tr(pri), should be understood in 
this way. 

In contrast, by a global measurement of observable ^l, we mean a measurement on the whole ensemble: 
we fetch every molecule in the ensemble and measure its fl. The result of the global measurement of Q will be 
the sum of all these individual measurements. The fluctuation of the global measurement of is understood in 
the same way as that for a quantum system: one prepares many copies of the same ensemble, he makes global 
measurement of il on each of these ensembles. The expectation value of the global measurement will be the 
averaged value of these global measurements, and the fluctuation is the standard deviation from this averaged 
(or mean) global measurement. We denote the global measurement by a subscript E. The expectation value of 
a global measurement of fl is denoted by {ft)E- For instance, for an ensemble of nuclear spins, the z-component 
of the spin operator for the whole ensemble is the sum of all the individual nuclear spins, 

N 



One can obtain its value by measuring the total magnetic polarization of the whole ensemble, thus sometimes 
one can get the global measurement result without measuring each of the individual molecule. 

The expectation value of a global measurement and the sampling measurement is related by the following 

{n)E=^N{n), (1) 

where N is the number of molecules in the ensemble. 

I . Consider the case in 



Von Neumann defined a density matrix and mixed and pure states in the following way 
which one does not know what state is actually present in a system, for example when several states |(/>i), |</'2), 
• • • with respective probabilities wi, W2, ■ • •, (wi > 0, W2 > 0, • • •, wi + W2 + • • • = 1) constitute the description 



of the system. Then the expectation value of a physical observable in the system is 

(n) = J2 w„{^nM(j)n) = Tiipn). (2) 

n 

The density operator, or density matrix 

P ^^Wn\(l)n){(t)n\, (3) 
n 

is used to describe the state of the quantum system under study. It characterizes the mixture of states just 
described completely with respect to its statistical properties. 

When only one Wi is equal to one and the rest Wi^s are zero, then the quantum system is described by a single 
wave function, and the system is said to be in a pure state. When more than one Wi's are nonzero, the state is 
called a mixed state. 

The above notion is commonly accepted. However, in practice there are different scenarios where the density 
matrix can be applied. For example, the abstract definition of density matrix is applied to an ensemble with N 
molecules. If in an ensemble of N independent molecules, there are Ni molecules in state \ipi)^ then the ensemble 
is described by the following density matrix 

m 

p = ^w^\^lj^){^il (4) 
1=1 

where Wi — Ni/N and m is the number of possible single molecule wave functions in the ensemble. 

Here the density matrix represents the "state" of an averaged molecule from this ensemble: a molecule in this 
ensemble has Wi — Ni/N probability in state Wi is the probability a molecule is in state if one fetches 
a molecule from the ensemble randomly. If m > 1, the ensemble is said to be a mixed ensemble, or simply 
a mixed state. If m = 1, the ensemble is an ensemble in which all molecules are in the same state, and the 
ensemble is called a pure ensemble, or simply a pure state. 

This definition of mixed state is widely adopted in textbooks and monographs of quantum mechanics, for 
example in Shankary|, Merzbacher and so on. d'Espagnat has suggested to use the term, proper mixture 



for this type of object 



and this practice is adopted in this article. 



^ and PreskiU ^ 



For example, in one of Peres's 



This definition of mixed state is also used by Asher Peres 
example, there are two different ensembles that have the same density matrix p — ^l2, where I2 is a 2 by 2 unit 
matrix. One ensemble is prepared by preparing each qubit in either state |0) or |1) according to the result of a 



coin-tossing. The other ensemble is prepared by making the state of each qubit in either 

|±) = (|0)±|l))/2, (5) 

also according to the result of a coin-tossing. Though each individual qubit from both ensembles is in a definite 
quantum state, but on average, a qubit from the first ensemble has half probability in state |0) and |1) 
respectively, namely if one fetches a qubit from the prepared ensemble, there is half probability the qubit is in 
state |0) and another half probability in state Similar statement could also be made for the second ensemble. 
One should note that in Peres's example, the ensemble is defined to have infinite number of molecules. In this 
article, we talk of ensembles with a finite number of molecules. 

6. Landau proposed the density matrix in a different way j2i- In a paper on damping problems in quantum 
mechanics. Landau introduced the density matrix to study coupled quantum systems. Suppose there are two 
systems, A and B, coupled together. The wave function for the coupled AB-system is 

I^ab) =^c„,|V„)IO. (6) 

n,r 

The expectation value of an observable pertinent to the A-system in this state is 

{n^) = {^AB\^^\i^AB) = = Tr(p^f^^), (7) 

where 

= (V'nl^^^lV'm), (8) 

and 

PA.nra = ^ C*^Cm,. = Tr^ ( | V^As) | ) ■ (9) 

r 

The density matrix pA describes the "state" of the A-system. When the rank of the density matrix is greater 
than 1, the state is also a mixed state. This type of mixed state is different from the proper mixed state, 
although they have the same mathematical structure. d'Espagnat suggested to use the terminology, improper 
mixture, for this type of object. We also adopt this in this article. 

7. There are two different situations that a probabilistic description is required. In one case, the description is due 
to our lack of information about the quantum system, but the quantum system itself is actually in a definite 



quantum state described by a wave function. For instance in quantum key distribution 



, Alice prepares a 



sequence of qubits, each qubit is randomly prepared in one of the following states: |0), |1), |+), |— ). To anyone 
other than Alice, each photon is in a proper mixed state described hy p ~ ^^2- But to Alice, it is in a definite 
pure state because she has prepared it herself. 

In another situation, a probabilistic description is the best one can provide. Suppose that two qubits A and B 
are in a Bell-basis state. 



\^ab) = \I^{\OaOb) + \IaIb)}. (10) 

If one wants to describe the state of qubit A without mentioning qubit B, the best one can do is to use the 
density matrix = ^12. If one measures the state of qubit there is half probability each to obtain |0) or 
Anyone, without exception, will get this result. This is in contrast to the previous example in which Alice has 
full knowledge of each qubit. 

8. There is an essential difference between the proper mixed state and the improper mixed state. In the proper 
mixed state, each individual molecule in an ensemble is in a definite quantum state. In an improper mixed 
state, the molecule is actually in such a state with the proposed probability distribution. In the example when 
qubits A and B are in a Bell-basis state, = {\OaOb) + |1a1s))/\/2, the reduced density matrix for qubit 
A is = il2. In this state the qubit A is actually in state |0) and |1) simultaneously. 

9. The mathematical formalism of the density matrix are well-known|^. Its equation of motion in the Schoedinger 
picture obeys the Liouville equation 

^h^^[H,p]. (11) 
The density matrix p is positively definite. Its properties are summarized as follows 

Tr(p) = 1, (12) 

= p, (13) 

Tr(p2) < 1. (14) 

It should be pointed out that wave function description can also be equivalently expressed in terms of a density 
matrix, 

P=l*>(*l- (15) 

The density matrix in this case of a pure state is characterized by the equality Tr(p'^) = 1. The density matrix 
with Tr(p^) < 1 is a mixed state, either proper or improper. 
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3. THREE DIFFERENT DENSITY MATRICES AND MIXED STATES 



3.1 The FULL DENSITY MATRIX for a single QUANTUM SYSTEM 



10. The state of a single, closed quantum system is always described by a state vector \ip). This is dictated by 
the basic postulate of quantum mechanics 9j: "the state of the particle is represented by a vector \ip{t)) in a 
Hilbert space." The density matrix for such a system is p = \ip){'ip\. In this case, the density matrix description 
and the wave function description are strictly equivalent, and the unphysical overall phase in the wave function 
description is naturally eliminated in the density matrix description. It contains ALL the physical information 
about the system. We call the density matrix in this case a full density matrix. 

A single, closed quantum system may have complicated structures. For instance a Bose-Einstein condensate in 
an isolated environment is considered as a single closed quantum system. An isolated atom with many electrons 
and nucleons is another example. A quantum system of N 2-level molecules(or qubits) could also be viewed as 
a single quantum system if we want to describe all the details of the whole N qubits. The density matrix for 
such an A^-qubits system is of 2^ x 2^ dimension. 

Though the state of such a system is always pure, a mixed state description for such a system is possible when 
one considers only the average properties of the system over a period of time, or over some parameters. For 
instance if the system is under rapid change in time and only average property of the system over a relatively 
long period of time is concerned, the state of the system could be described by a density matrix obtained by 
averaging over the time 



where X{t) is some probability distribution function. If the resulting density matrix describes a mixed state, 
it is a proper mixed state because we have actually traded the space averaging in an ensemble with the time 
averaging for a single molecule over some time. 




tl 



(16) 



3.2 The COMPRESSED DENSITY MATRIX for an ENSEMBLE 



11. We call the single- molecule density matrix over an ensemble of N independent molecules as a compressed 
density matrix. If the compressed density matrix is pure, it indicates that all the single molecules in the 
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ensemble are in the same quantum state. If the compressed density matrix is a mixed one, it represents 
an ensemble with molecules in distinct quantum states. As stressed earlier, each individual molecule in the 
ensemble is in a definite quantum state, rather than in different states with respective probabilities. A mixed 
state arising from the average over an ensemble is called a proper mixture by d'Espa gna tQ as we mentioned 
earlier. Mixed state from this average is called the ignorance interpretation by Cohen |l4| . 



3.3 REDUCED DENSITY MATRIX in composite QUANTUM SYSTEM 



12. The density matrix obtained by tracing out partial degrees of freedom of a coupled system is called reduced 
density matrix. Suppose there are M constituent molecules in a coupled quantum system. The wave function 
of the coupled system is \'i'AiA2...AM)- The reduced density matrix for molecule Ai can be obtained by tracing 
out the degrees of other molecules 

PAi = TrAiA2...Ai_iA,+i...AM(l*AiA2...AM>(*AiA2...AM|)- (17) 

One can also define reduced density matrix for two or more molecules by tracing the degrees of freedom of the 
remaining molecules, for instance, the reduced density matrix for molecules Ai and Aj is defined as 

PAiAj = TrAiA2...Ai_iA,+i...Aj_iA, + i...A„(|*AiA2...A„)(^'AiA2...AMl)- (18) 

The reduced density matrix has also been defined for identical particle systems. For an M-identical-particle 
system, one can define one-particle, two-particle, ... reduced density matrices as follows ^ 

= Tr,,,3...,^,(p), 
= Tr,3...,„(p), 



P^''-'^ = Tr,,,(p), (19) 



where p is the density matrix for the coupled M-particle system. The whole M-particle system is treated as a 
single quantum system. The subscript ik is the single-particle state index. These definitions are useful in the 
study of quantum entanglement in identical particle systems [ 

It is worth noting that an M-particle coupled quantum system is different from an ensemble of M independent 
molecules. Usually the particles in an Ai^-particle coupled quantum system are interacting and correlated. 
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However the molecules in an ensemble described by a compressed density matrix are independent and non- 
interacting. The coupled M particle system should be described by a Af-particle density matrix, while an 
ensemble of M molecules is described by a single molecule compressed density matrix. For example, a M-qubit 
system is described by a 2^^ x 2^^ density matrix, and an ensemble of M qubits, such as M spin-1/2 molecules is 
described by a 2 x 2 density matrix. In an NMR ensemble of 7-qubit molecules, there are usually N = 0(10^^) 
molecules. However the compressed density matrix for the ensemble is 2^ x 2^, rather than 2^^ x 2^^ 18l |. 

An improper mixed state usually could not be identified with a proper mixed state with the same density 
matrix|5j, as we see they represent different objects though their mathematical descriptions are the same. This 
type of mixed state was interpreted by Cohen as the ancilla interpretation jljj . 

13. It is useful to make a comparison between a single qubit and a coupled 2-qubit system to demonstrate the 
difference between a pure quantum state and an improper mixed state. Suppose a single qubit A is in state 

\+^a)^^{\0a) + \1a)}. (20) 

Consider another system in which qubit A and qubit B form an entangled state 

1 



AS/ - Y 2{|0a0b) + |U1b)}. (21) 

If one measures az on qubit A, he will obtain \0a) and \ Ia) with 50% probability in each case. But if he measures 
ax on qubit A, he will obtain with certainty | + xa) in the first case as it is an measurement on an eigenstate, 
whereas he will still have 50% probability to obtain | + and \ ~ xa) respectively in the latter case. 



4. THREE SIMPLE EXAMPLES 



14. We study three ensembles of qubits, that will be all described by the same mathematical density matrix p = il2 
when N approaches infinity, to illustrate the different mixed states and their relationship. Suppose there are N 
pairs of qubits in state given by Eq. (|21|) . All the A qubits are given to Alice and all B qubits to Bob. They form 
ensembles Ai and Bi respectively. Hence each A qubit molecule in Ai is in an improper mixed state described 
by reduced density matrix pA = \l2: and similarly each B qubit molecule in ensemble Bi is also in an improper 
mixed state with the same reduced density matrix. It should be pointed out that all the A qubit molecules 
in ensemble Ai are in the same state, there is no difference between the states of the different molecules in 
ensemble Ai. The same conclusion can also be said about ensemble Bi. If one looks at the compositions of the 
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ensembles, one could say that the ensembles are "pure", as they are made of molecules all in the same "state", 
though the "state" is not a pure quantum state described by a single state vector. This is a generalization of 
the usual pure state in which all the molecules are in the same quantum state described by a wave function. 

In the second example, one performs cfz measurement on each qubit in ensemble Bi. Because of the collapse 
of state in quantum mechanics, each qubit in B\ collapses into either |0) or |1), with equal probabilities. 
Consequently, the corresponding A molecules in ensemble Ai also collapses into state |1) or |0) respectively. 
The ensembles after the measurement are called A2 and B2 respectively. In ensemble B2, each qubit is in a 
definite quantum state and the preparer of the ensemble who performs the az measurement knows all these 
information. He also knows the state of A molecules in ensemble A2 as he can infer from the prior entanglement 
between qubit A and B and the result of his measurement of qubit B. Suppose after the measurement, the 
number of qubits in state |0)(or |1)) is N/2 — Nsipr: N/2 + Ns), then the compressed density matrix of A2 is 



/ 1 



PA2 



2 Af " 



\ i + 

When N approaches infinity, the compressed density matrices for qubit A2 and B2 all become p = |l2. Ng is a. 

random number, can be positive and negative or zero, and is proportional to -s/lv according to statistics. 

In the third example, one performs measurement on each qubit in ensemble Bi. Again due to the collapse 
of state in quantum mechanics, each qubit in Bi collapses into either \+ x) or \ — x), with equal probabilities. 
Consequently, the corresponding B molecules in ensemble Ai also collapses into state \ —x) or \+x) respectively. 
The ensembles after the measurement are called A3 and B3 respectively. In ensemble B3, each qubit is also in 
a definite quantum state and the preparer of the ensembles who has performed the ax measurement knows all 
these information. He also knows the state of A molecules in ensemble A3. Suppose after the measurement, the 
number of qubit in \ + x) (or | — a;)) is N/2 — N^ (or N/2 + N^), then the compressed density matrix of A3 is 



(22) 



/ 1 _N^\ 
2 N 

\ N 2 / 



(23) 



Again the compressed density matrices for qubit A3 and B3 are all approaching p = il2 when A'' approaches 
infinity. Here N^ is also random, could be positive, negative or zero, and is proportional to y/N. 

We stress here that when the number of qubit A^ is finite, the three ensembles are different and they are 
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distinguishable physically. When a measurement is performed on an ensemble, the compressed density matrix 
of the resulting ensemble is usually uncontrollable, and this is especially so when the number N is small. Even 
if one performs the measurement on the same observable, the resulting compressed density matrix is different at 
different times. Suppose one has N EPR pairs in state in Eq. H21|) . As usual, the A particles are given to Alice 
and the B particles given to Bob. When Bob performs Uz measurement on each of his B particles, he will get 
an ensemble with a compressed density matrix pg. Suppose he repeats the process a second time starting from 
state in Ea. (|21|l . he will obtain another ensemble with a compressed density matrix p'g. Usually ^ as 
can be seen in Ea. (|22|l . An extreme case is an ensemble Ai with N = 2. Before any measurement, the state of 
any of the two particles is an improper mixed state with density matrix, p = I/2I2. Suppose az is measured on 
each of the two particles. Then four different results may occur: 1) both qubits have spins up, and the resulting 
compressed density matrix is a pure state, 



1 


2) both qubits have spins down, and this gives a compressed density matrix, 

/ 



(24) 



(25) 



V 




1 

3) one qubit has spin up and the other spin down, and vice versa. In both cases, the compressed density matrices 
of the resulting ensemble are 

/ 1 



2 



(26) 



This third case is more probable, and there is 50% probability to occur, whereas case 1 and 2 each has 25% 
probability to occur. 

Though the three ensembles are distinguishable, they can not be exploited to realize faster-than-light quantum 
communication because one can not ascribe a unique compressed density matrix to a measurement because of 
the fluctuation in the cjuantum measurement. Impossibility of faster-than-light communication docs not imply 
the indistinguishability of ensembles having the same compressed density matrix, which has been discussed 
recently. 



jimg 
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5. DISTINCTION OF ENSEMBLES HAVING THE SAME COMPRESSED DENSITY MATRIX 

5.1 Present STATUS 

15. Penrose doubts the completeness of the density matrix description. In an example where A (I myself) on Earth 
and B (my colleague) on the Moon each hold a particle from an EPR pair, Penrose states that "However, it is 
hard to take the position that the density matrix describes reality. The trouble is that I do not know that I 
might not, sometimes later, get a message from the Moon, telling me that my colleague actually measured the 
state and found the answer to be such and such. Then, I know what my particles's state must actually be. The 
density matrix did not tell me everything about the state of my particle. For that, I really need to know the 
actual state of the combined pair. So, the density matrix is sort of provisional description, and that is why it is 
sometimes called FAPP(i.e., for all practical purposes)." 

Here Penrose refers to what we call the reduced density matrix for a composite quantum system of two constituent 
particles. In our later discussion, we mainly focus our attention to the compressed density matrix. 

16. It is commonly believed that ensembles with the same compressed density matrix are physically identical. We 
believe that this misbelief can largely be eliminated once a clear distinction between the compressed and the 
full density matrices is established. 

A complexity arises due to the indistinction of proper mixed state and improper mixed state. Here we are 
talking about the distinction of ensembles with the same compressed density matrix rather than different com- 
posite systems with the same reduced density matrix. At present, there is intensive research work on quantum 
entanglement based on reduced density matrix for composite quantum systems. The distinction of ensembles 
having the same compressed density matrix we address here is not concerned with coupled quantum systems. 

After clearing up those confusions in the terminology, there are still some authors who believe that ensembles 
with the same compressed density matrix are physically identical. 

17. Asher Peres stresses that ensembles with infinite number of particles having the same compressed density matrix 
could not be distinguished physically. He uses an example to illustrate this idea. In one ensemble with infinite 
number of photons, each photon is prepared in either the \H) = |0) state or \V) = |1) state according to the 
result of coin-tossing. Hence the compressed density of this ensemble is p = ^^2- In another ensemble of infinite 
number of photons, each photon is prepared in the left-handed circular or right-hand circular state according 
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to the result of coin-tossing, and the compressed density matrix is also p = ^12. Peres argued that there are 
no physical means to distinguish these two ensembles. He has elevated this into a fundamental postulate in 
quantum mechanics 

The (compressed ) p matrix specifies all the physical properties of a quantum ensemble (with infinite 
number of molecules). 

Here the words inside the parenthesis are added by the authors to conform to the terminology used in this 
article. 

At first sight, it may seem that this fundamental postulate excludes any possibility of distinguishing ensembles 
with the same compressed density matrix. However, it is not so because Peres' definition of ensemble is different 
from ours. Peres' ensemble contains an infinite number of molecules and it only exists conceptually as stressed by 
Peres y. Thus it does not exclude the possibility of distinguishing ensembles with a finite number of molecules 
having the same compressed density matrix. 

In fact, Peres stated that finite particle number ensembles with the same compressed density matrix can be 



physically distinguished 



3 



|. In this example, Peres constructed two ensembles with 1 million qubits each. In 
one case, half of the qubits are polarized in |0) state and the other half in state |1). In another case, half of 
the qubits are in state -^(|0) + i|l)) and the other half in state -^(|0) ~ With high probability, people 

can distinguish the two ensembles (assemblies in Peres's term). The probability of failure is proportional to 
where N is the number of molecules in the ensemble. Apparently, as N increases, the probability of failure is 
approaching zero. 

18. Preskill stresses that differently prepared ensembles with the same compressed density matrix are not distin- 
guishable on one hand, but he also points out different preparations of ensembles with the same compressed 



y. 



It is noted that Preskill's 



density matrix could be distinguished physically with some additional external help 
ensemble has a fixed number of molecules. Explicitly, Preskill proposed a method to distinguish the ensembles 
A2 and ^3 in section 4. He presumably assumes that the compressed density matrix of A2 and are the 
same, which is usually not true when the number of particles is a fixed number. Suppose Bob starts from 
ensemble Bi and makes a Ux measurement on his B qubit molecules and hence prepares the ensemble and 
A3. Alice can find out what ensemble her ensemble is with some help from Bob: Bob tells Alice the result of 
the measurement for each B qubit molecule, whether up or down, but he does not tell her what measuring basis 



15 

he has used. Ahce can choose any measuring basis, whether cr^: or cr^, to measure each of her A qubit. If Ahce 

chooses ax, her resuh will have perfect accordance with Bob's. If she chooses CTz, only 50% of the results agree 

n 

with Bob's. Hence Alice "does have a way to distinguish Bob's two preparation methods" (j] . 

19. d'Espagnat pointed out that ensembles having the same compressed density matrix can be distinguished phys- 
ically by observing fluctuations of some observables. He has illustrated this in an ensemble of N qubits 5|. In 
ensemble 5*1, a half of qubits are in state |0) and another half are in state |1) where |0) and |1) are the eigenstates 
of operator ci^, the Pauli matrix. In ensemble 5*2, a half of qubits are in state | + x) and another half are in 
states I — x) where | -I- x) and | — x) are the eigenstates of operator cTx- In both ensembles, the compressed 
density matrix is p = ^12. Although the average value of any observable f2, in our terminology in a sampling 
measurement, for both ensembles are all expressed as 

{n) = Tr(pf]), (27) 

the fluctuation of an observable will be different for ensembles with different compositions. Translated into our 
language, d'Espagnat actually refers to the fluctuation of a global measurement is different for different ensemble 
compositions with the same compressed density matrix. d'Espagnat used the following observable 

N 



=^a.(i), (28) 



where the summation over i is made over all the molecules in the ensemble. For ensemble ^i, the fluctuation of 
Ez is zero, 

(AS,)i5 = 0, (29) 

because each molecule is in an eigenstate of az and every measurement of on the whole ensemble will give 
identical result. However for ensemble 5*2, the fluctuation is 



(AE,)i; = VA^, (30) 

namely if one measures Sz on the whole ensemble, different result may be obtained at different times. By 
observing the ffuctuation one can distinguish between the two ensembles with the same compressed density 
matrix. 
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5.2 ENSEMBLES with the same FULL DENSITY MATRIX are identical 

20. Suppose one treats the whole ensemble having N molecules as a single quantum system, the ensemble is described 
by a full density matrix. If each molecule in the ensemble is independent and non-interacting, then the state of 
the ensemble is given by 

IVb) = |Vi(l))|V'i(2)> • • • |V'i(^i)>l^2(iVi + 1)) • • • \MNi + iV2)) • • • , (31) 

where the number inside the round bracket is the labelling of the molecule and the subscript refers to the state in 
which the molecule is in. The state {iPe) is the product of states of individual molecules. The density matrix 
for the whole ensemble of iV molecules is p = |V'b)(V'b|- Since the wave function contains all the information 
about the system, the corresponding density matrix completely describes the state of the whole ensemble. We 
conclude that ensembles with the same full density matrix are physically identical. In this case, the density 
matrix is always one representing a pure state of a quantum system with N molecules. If the full density 
matrices of the two ensembles are the same, then the two ensembles are physically identical. 

5.3 ENSEMBLES having the same COMPRESSED DENSITY MATRIX but with different PARTICLE 

NUMBERS are physically distinguishable 

21. First, let's discuss the case of two ensembles having the same compressed density matrix but with different 
molecule numbers. With ideal measuring sensitivity, one can measure up to the accuracy of single molecule. 
Then from the expectation value of a global measurement alone, one can distinguish the two ensembles. The 
expectation of any observable with the form 

N 

Q = "^n{i}, (32) 

will be sufficient to distinguish the two ensembles, because the expectation of observable O is 

= 7VTr(pO). (33) 

Because different molecule number N will give different expectation value, one can distinguish these ensembles 
from the expectation value directly. This is easy to understand if we compare two different amount of samples 
in nuclear magnetic resonance experiment. In NMR, the measured quantity is the total magnetization, and is 
equivalent to the signals in coil. The larger the sample, the stronger the measured signal. 
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5.4 General EXPRESSION for the FLUCTUATION of OBSERVABLES in GLOBAL MEASUREMENT 

22. Generalizing the idea of d'Espagnat, ensembles having the same compressed density matrix could be distin- 
guished by inspecting the fluctuation of some observables in global measurement. By this method, we can 
distinguish ensembles upto different physical composition: the number of molecules Ni in state {ipi)- 

Suppose there are iV molecules in an ensemble with Nk molecules in state {ipk), and fc = 1, • • • , m. The fluctuation 
of the observable 

N 

= (34) 

i=l 

where i is the molecule index in the ensemble, is 



{An)E = ^NiTTin^p)-Y,Nk{M^\^pk)^ (35) 

where the summation is over the possible single molecule states in the ensemble. It is apparent that the second 
term under the square-root is composition sensitive: ensembles with different compositions will give different 
results. This can be comprehended in the following way. For a single molecule in state iV'fe): the fluctuation of 
n squared is 

iAn)l = (v-fcirj'l^fc) - {i^k\n\iJkf, (36) 

and contribution from Nk such molecules is Nk{An)1, and the total contribution from all the molecules in the 
ensemble is J2k ^fc(Arj)|. With some simple calculation, this gives the result in Eq. H35II . 

This result can be derived directly if one recalls that the state of the whole ensemble is Eq. (|31|l The fluctuation 
of observable fl, namely the fluctuation of the global measurement is 



iAn)E = ^J{^E\n^~m))l\^|;E). (37) 
The expression inside the square-root can be calculated and it is 

{^peI (5]f7(z)2+^f](z)f7(jvE(^^i"«i^^^)'-E(^^i"(W)i^^)| 1^^) 
fc fc 

where i runs over the molecules in the ensemble and k runs over all possible single molecule states in the 
ensemble. 
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5.5 Composite QUANTUM SYSTEMS with the same REDUCED DENSITY MATRIX 

23. Different composite quantum systems may produce tfie same reduced density matrix. They could be distin- 
guished if one measures quantities involving the whole composite quantum system. For instance let us consider 
two systems of an EPR pair A and B, one in state 

\^\b)=^\IIwaOb) + \IaIb)}, (39) 

and another in state 

\^Ib)^\IIwaOb)-\IaIb)}. (40) 

In both cases, the reduced density matrix is p^i = -^12. By measuring the observables of molecule A alone, 
one can not distinguish these two different composite quantum systems. However, the two systems will be 
distinguished the moment one makes a joint Bell-basis measurement. 

6. THE NATURE OF ENSEMBLE NUCLEAR MAGNETIC RESONANCE QUANTUM 

COMPUTING 

24. Ensembles with the same compressed density matrix can be distinguished physically is very important in quan- 
tum information processing. One important issue is the nature of room-temperature NMR quantum computing. 
It has been widely accepted that NMR quantum computing is not quantum mechanical because there is no 



entanglement in each step of the NMR quantum computing 



22| . We point out here that this conclusion is based 



solely on the unjustified belief that ensembles having the same compressed density matrix are identical. 

The major argument in Ref. 2^ is to rewrite the effective density matrix, which is the compressed density matrix 
of the effective pure state obtained after the effective pure state technique, 

'1 - 



Id + epeff, (41) 

where Pe// is the effective pure state density matrix, and e is a number about 10^^ large and Id is d x d unit 
matrix, into a convex decomposition of product-state density matrices, 

p^Y.^4p^'s>---p^^y ' (42) 



19 



Pi P2 P3 P4 P5 Pe 
Pi e -e 
P2 -e e 
P3 e e 

P4 -e e 
P5 e -e 

Pe -e e 



TABLE I: Cij for the effective pure Bell-state. A space means that coefficient is zero. 

where Ci is the decomposition coefficient and {pi • • • Pa/)* is the density matrix representing the product of 
M qubits. This convex decomposition is possible when e is small and the qubit number M is small. Since the 
product state does not have entanglement, then it was concluded by the equivalence of ensembles having the same 
compressed density matrix that the current NMR quantum computing does not have quantum entanglement, 
and hence is not genuine quantum mechanical. 

From the previous discussions in this article, we see that the conclusion is based on an unjustified belief. 
Though the compressed density matrix of the effective pure state could be rewritten into convex decomposition 
of products state density matrices, they are not physically equivalent. As a simple example, we show that the 
effective Bell-state ensemble and its product state decomposition in Ref.|2^ can be distinguished physically. 
The effective Bell-state compressed density matrix is 

P =^—^h + tPBelU (43) 

where pij^u = (|00) + |11))((00| + (ll|)/2. The effective Bell-state describes an ensemble with eN molecules in 
Bell-state (|00) + |ll))A/2 and (1 - e)N/A molecules in each of the calculating-basis states |00), |01), |10) and 
1 11). According to Ref.l2^, it can be decomposed into 

P = T.\(l + ^^i)P^^Pi' (44) 

where Pi — (I + <7i)/2 for i = 1,2,3 which represents a pure state polarized along the three axis x, y and z 
respectively, and P,; = (1 — CTi)/2 for i — 4, 5, 6 which represents a pure state anti-polarized along the three axis 
X, y and z respectively. The coefficient dj are given in Tabled The two different compositions of the effective 
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Bell-state can be distinguished by observing the fluctuation of the operator 

Sz. =5I'^i-(*)'^2.(j)- (45) 

i 

For the ensemble implied by (|43|l . the fluctuation is {/^Tizz)i.E — ^VN, whereas for the product state expansion 
gH), the fluctuation is (AS^^)2,_e = 

Hence the two ensembles are distinguished though they all have the same compressed density matrix. The de- 
composition of the effective Bell-state compressed density matrix into a product state compressed density matrix 
can not be used to infer the conclusion that no entanglement exists in the effective Bell-state. Entanglement is 
the property of the individual molecules, not the ensemble as a whole. As in Ref. 2^, even if all the molecules 
in an ensemble are in the same quantum state, there is no entanglement between different molecules. Hence the 
effective pure state should be considered as what it is physically, and there are indeed quantum entanglement 
within each molecule, though only a small portion of the molecule from the effective pure state contribute to 
the final read-out in the NMR detection signal. 

25. In addition, there is an unnaturalness difficulty in the product state decomposition description of the NMR 
quantum computing used in Ref. [2^. In the effective pure state NMR quantum computing, the quantum 
computation is performed on the eN molecules in the effective pure state. In the product state decomposition, 
if only single qubit operations are performed one can still retain a simple picture. We can visualize that there 
are several "tubes" in the NMR sample. For instance in a two-qubit system case, there are 36 "tubes" in the 
sample, each tube contains some number of molecules in a product state Pi ® Pj. A single qubit operation just 
changes the single molecule state into another state so that one can still use the same "tubes" to describe the 
state of the ensemble. However once an entangling operation such as the controUed-NOT gate is performed, the 
number of molecules within each composition, that is in a "tube", will change. For example, if we perform a 
CNOT gate on effective Bell-state (I43f) . this changes the ensemble into an effective | -I- x, 0)-state, that is, the 
first qubit is in | + x) state and the second qubit is in |0) state, 

= ■^^^/4 + e| +a;,0)(+x,0|, (46) 



which is natural and clear. However in the product state expansion of Ref. 22], this density matrix becomes 
where the coefficients C[ ^ are given in Table |HJ We see that there is a reshuffling of the number of molecules 
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TABLE II: C'ij for the effective pure state | + a;,0). A space means that coefficient is zero. 

in different states. In order to get the same averaged result of the effective | + a;, 0) state, one has to reshuffle 
the distribution of the molecule numbers in different "tubes". While for non-entangling operations, there is no 
need for such a reshuffling because the single-qubit operation can be performed on each composition of states 
directly, and their effect is to change state Pi to state P/ = UPiW. 

One objection to Eq. (|36|l is that each molecule in an ensemble is in a mixed state itself, even at a given instant. 
In this case, the fluctuation of observable is 



A^E = V^Tr(p02) - iV(Tr(pO)2. (48) 

In this case, all molecules in the ensemble are equivalent and no physical means can distinguish them. In fact, 
there is no difference at all, because the compositions of the ensembles are the same. This has been used by 
some authors to criticize d'Espagnat To this objection, we stress that it is impossible to prepare such an 
ensemble using classical probabilities. For instance, in the example given by Peres the state of a photon is 
prepared according to the result of a coin-tossing, either along z or along — z. It is true that each photon has 
50% probability to be prepared in state |0) or |1). But each photon is prepared in a definite state, not in a 
mixed state. This is exactly the case in the BB84 quantum key distribution protocol ^ where the state of a 
qubit has 25% probability in each of the four states \z), \ — z), \x) and | — a;)- But it is certainly in one of the 
four possible states. 

One possible way to realize such a scenario is that each molecule is in an improper mixed state. For instance in 
an ensemble of molecules, each molecule contains two qubits A and B, and the two qubits are in an entangled 



icitly 



22 

state, say 

|^ab) = ^ (100)^5 + |11)as). (49) 

If one looks at the state of the first qubit in each molecule in the ensemble, the resulting state is an improper 
mixed state. Then in this case, by observing the fluctuations of observables for qubit A alone one can not 
distinguish another ensemble with the same compressed density matrix where A and B are in state 

IVW = ^(|01)ab + |10)ab). (50) 

n„_, . „o. . .M. c»,»Ma „a c.„». .p,.p 

explained that in effective pure state there are some numbers of molecules in a definite quantum state[] 
This is a valid approximation for NMR quantum computing because nuclear spins are well isolated from the 
environment, hence it has a long decoherence time. For those nuclear spins that are not used in quantum 
computing, decoupling pulses are usually used to disentangle them from the working nuclear spins so that they 
can not form entangled state that will make the working qubits in an improper mixed state. Thus a valid picture 
of the ensemble used in NMR quantum computing is that the ensemble contains lots of molecules, there are 
certain number of molecules in each of the quantum states of the molecule. The spin-relaxation time Ti can be 
viewed as the time period a molecule remains at a given quantum state. This is supported by he Gorter formula 

y 

1 _ 1 En.m WnM^rn - E^f 

where Wm,n is the transition probability rate from level m to level n. Hence Ti is the energy weighted time a 
molecule remain at a given quantum state. 

What is then the significance of studying quantum entanglement directly from the density matrix? The study of 
entanglement for the compressed density matrix directly is of significance because for those compressed density 
matrices that can not be decomposed into product states, their microscopic entanglement are already reflected 
in the expectation value. However for ensembles, it is more important to study the entanglement of the pure 
state of single molecules because it is the place where quantum entanglement exists. As we see, even though 
the expectation value does not show the entanglement, one can still exploit this precious resource through some 
method. On the other hand, for improper mixed state, the study of entanglement based on density matrix is of 
importance. 
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7. SUPPORTING EXAMPLE FROM QUANTUM COMMUNICATION 



27. Though fluctuations of observables in ensembles with infinite number of photons can not distinguish different 
preparations, there are stiU methods to distinguish them if one has individual access to the photons in the 
ensemble. The quantum key distribution is a good example. First we briefly describe the BB84 quantum 
key distribution scheme}]^. In the BB84 scheme, Alice sends a sequence of photons to Bob. Each photon is 
randomly in one of the four states: |0), | + x), | — a;) respectively. Hence the compressed density matrix of 
the ensemble, with inflnite number of photons theoretically, is p = ^l2- Then Bob chooses randomly one of the 
two measuring devices, the cr^-basis or (Tjj-basis to measure the state of the photons. After the transmission, 
Alice and Bob publicly compare their measuring-basis used for each photon, and they retain those events in 
which they choose the same measuring-basis. From these retained events, they choose a sufficiently large subset 
of events and publicly compare the measuring results of these events. If the error rate is high, they conclude 
that there is Eve during the transmission and discard the results. Otherwise they proceed to the post-processing 
stage and get a final secret key. 

If the compressed density matrix had described all the physical properties of the ensemble, we would run 
into difficulty in the security of quantum key distribution. The same compressed density matrix could also be 
produced by Alice with only two states |0) and |1): she randomly prepares the photons in one of these two states. 
In this way, the compressed density matrix is the same p = \^2- Immediately we know that this replacement is 
invalid: the simplified scheme is unsafe, for if Eve uses the Cz measuring device to measure every photon, she 
can steal every bit of information. Then she can resends a photon with the same state to Bob. In this way, her 
action could not at all be detected. Hence the scheme becomes insecure. This conclusion would not alter if we 
had used infinite number of photons in the ensemble. 



8. SUMMARY 



28. We have examined in details the definition of ensembles in quantum mechanics and the different definitions of 
mixed states. We have proposed to use different terminologies for the density matrix in different situations: the 
full density matrix, the compressed density matrix and the reduced density matrix. We have explicitly shown 
that ensembles having the same compressed density matrix but different compositions can be distinguished 
physically by observing fluctuations of observables in the whole ensemble. With this conclusion at our hand. 
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we studied the problem of the nature of NMR quantum computing. We have pointed out that the conclusion 
that there is no entanglement in the current NMR quantum computing experiment is based on an unjustified 
belief that ensembles having the same compressed density matrix are physically equivalent. This conclusion is 
also supported by the security aspect of quantum key distribution where the density matrix p could not specify 
all the physical properties of an ensemble, even if the number of molecules in the ensemble is infinite. A brief 
summary is given in Table UlTI 
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TABLE III: A summary 



Pull density matrix 


Object described 
Matrix Name 
Object Name 

Remark 


A single quantum system, but may contain many molecules 
tull density matrix p/ = 

A pure state. The object is always described by a single wave function, or equivalently 
a full density matrix for a pure state 

The description is complete. Different preparations with the same full density matrix 
are physically identical. 


Compressed density matrix 


Object described 
Matrix Name 

Object Name 
Remarks 


An ensemble of N independent molecules with Ni molecules in state 
Compressed density matrix pc = jf-\il)i)(ipi\. pc describes the state of an averaged 
molecule from this ensemble. The full density of the whole ensemble is p/ = pi iS> 
P2 iS> • • • <S> pjv, where the state of different molecules can be the same. Though two 
molecules may be in the same quantum state, they are distinguishable, for instance 
by their positions in space. The relation between the full density matrix and the 
compressed density matrix is Pc = jf'^iPi, where pi is the full density matrix for 
molecule i in the ensemble. 

The object described by pc is called a proper mixture. 

One can not tell the difference between two ensembles having the same compressed 
density matrix by sampling measurement or by the average value, but they can be 
distinguished by observing the fluctuations of observables regarding the whole ensem- 
ble. 


Reduced density matrix 


Object described 
Matrix Name 

Object Name 
Remark 


Composite quantum systems 

Reduced density matrix which is obtained by taking trace over other degrees of free- 
dom from the full density matrix of the whole composite system, 
improper mixture. 

By observing the observables of only a part from the composite system, it is impossible 
to distinguish systems with the same reduced density matrix. But the composite 
systems can be distinguished if the whole composite systems are measured. 



